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Abstract. A method for computing global minima of real mul- 
OO ' tivariate polynomials based on semidefinite programming was de- 

veloped by N. Z. Shor, J. B. Lasserre and P. A. Parrilo. The aim 
of this article is to extend a variant of their method to noncom- 
r ) , mutative symmetric polynomials in variables X and Y satisfying 

^ ■ YX - XY ^ 1 and X* ^ X, Y* = -Y. Global minima of such 

, polynomials are defined and showed to be equal to minima of the 

• spectra of the corresponding differential operators. We also discuss 

d , how to exploit sparsity and symmetry. Several numerical experi- 

ments are included. The last section explains how our theory fits 
into the framework of noncommutative real algebraic geometry. 

> 

1. Motivation 

(N 

CN I One of the most popular methods for computing global infima of real 

^ ■ multivariate polynomials is the method of sums of squares relaxations. 

The idea goes back to N. Z. Shor (his original papers are summarized 
O ■ ill [Ml Chapter 9]) and it was further developed by J. B. Lasserre [13] 

and P. A. Parrilo [21]. See [H] for an extensive survey and [2l| [8] 
for implementations. In this section we will present a variant of this 
method and in the next section we will extend it from polynomials to 
5^ ! polynomial differential operators. Later sections are concerned with 

improvements of the basic method and numerical experiments. 
For a given polynomial / G IR[X], X = (Xi, . . . , Xd), write 

(1) inf / := inf {/(x) | x G M'^} = sup{/i G M | / - > on M*^} 

for its unconstrained global infimum. By convention inf / = — oo if / is 
not bounded from below. Note that the polynomial Xf -|- (1 — XxX-^f 
is bounded from below but it does not attain its infimum. 
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Write Yl for the set of all sums of squares of polynomials from 

Clearly, every element from nonnegative on W^. In 

particular this is true for the polynomial s := 1 + Xf + . . . + X^. On 
the other hand, the Motzkin polynomial 1 + XfX^ + XfX^ - SXfX^ 
is nonnegative on W^, d > 2, but it does not belong to 

Let us consider the following sequence of approximations of inf /: 

(2) := sup{/x G R I s\f - ^) G m[X]'}, 

where k e N := {0,1,2,...}. Since Xl^^l:^]^ is closed in the finest 
locally convex topology of ]R[X] by j29^ Corollary 11.6.4], we can replace 
sup by max in ([2]). 

Proposition 1. If f E satisfies the following assumption: 

/*) fn{x) > for every nonzero x eM.'^ where n = 2m is the total 
degree of f and fn is the n-th homogeneous part of f . 

then lim fik{f) = inf /. 

Remark. The assumption (*) is sufficent for the existence of global 
infimum but it is not necessary. (A necessary condition is that n = 2m 
and fn{x) > for every nonzero x G M.'^.) 

Proof. It is clear from ([1]) and ([2]) that inf / > ^k{f) for every k. Since 
s^{f - /i) G implies s^+\f - /i) G we have by P 

that > AiA;(/) for every k. Finally, for every e > 0, / — inf / + £ 

is strictly positive on W^. By Theorem [1] we can find /c^ G N such that 
s^^{f-m{f + e)eY.nXf- Hence, /ifc,(/)> inf/ -£ by ©. □ 

The following dehomogenized version of a theorem of Reznick [26] 
was used in the proof, see the comments after Theorem 5.5.2 in [T7] . 

Theorem 1 (Reznick 1995). // / G ffi[X] satisfies f{x) > for every 
a; G M'^ and the property (*) then s^f G ^M[X]^ for some A; G N. 

Finally, we would like to convince the reader that the numbers fJ^kif) 
can be effectively computed. We will do so by reformulating the defini- 
tion of fJ'kif) as a semidefinite program, i.e. an optimization problem: 

(3) minimize tr(CZ) subject to tr(AjZ) = 6j and Z > 

where C, Ai are given real symmetric matrices of the same size, bi are 
given real numbers and Z is an unknown real symmetric matrix of the 
same size as C, Ai. 

Let Vk be a vector of all monomials of total degree less or equal to 
i deg(sV) = k+m. Its size is = (''+7+'^) . Every element of ^ M[X]2 
of degree < 2{k + m) (in particular s^{f — /x) for every /i < fik{f)) can 
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be written in the form v^Zvk where Z is a positive semidefinite real 
symmetric matrix of size n^. Therefore, by 02]), we can express Hk{f) 
as the solution of the following optimization problem: 

, . Find fikif) = niaxjr/i where = {{ft, Z) | /i G M, 
Z G M„,(M), Z^ = Z, Z>0, s'if -fi)= vlZvk). 

Note that (|4]) is not exactly a semidefinite program as defined by ([3]) 
but it can easily be converted into one by: 

• eliminating /i using the linear relation obtained by comparing 
constant terms in s^{f — fi) = v'^Zvk, (we get /o — /i = Zn if 
/o is the constant term of / and the first component of Vk is 1), 

• writing the linear relations between Zij obtained by comparing 
coefficients in s''{f — fo + Zu) = vJZvk in the form tr(AjZ) = bi, 

• replacing /ifc(/) = max/x by /o-/ifc(/) = minZn = mintr(CZ) 
where Cu = 1 and other Cij are 0. 

Remark. Recall that the Newton polytope N{f) of a polynomial / = 
^CqX" G M.[X] is the convex hull of its support {« | Cq, 7^ 0}. The 
main property of Newton polytopes is N{fg) = N{f) + N[g). (It is 
proved by showing that both sets have the same extreme points, i.e. 
extiV(/^7) = extiV(/)+extiV(^) = ext{N{f) + N{g)).) If / = E^?^, the 
property implies that N{gi) C ^N{f) for every i. Hence, the vectors 
Vk from (jlj) need not contain all monomials of degree < |degs'^/ but 
only the monomials from ^N{s^{f — fi)). 

2. Polynomial differential operators 

Our aim is to develop a similar theory for hermitian elements of the 
d-th Weyl algebra W((i). Recall that W((i) is the unital complex *- 
algebra with generators X^, Yi, defining relations V/X^ — X^Yi = 5u 
and involution XI = X^, Yj* = —Yi, where k,l = l,...d. We will 
write W{d)h for the set of all elements u G yV{d) such that u* = u and 

VV((i)^ for the set of all finite sums of elements u*u where u G yV{d). 

The Schrodinger representation ttq is the ^-representation of yV{d) 
acting on the Schwartz space 5(M'^) considered as dense domain of 
L^(R'^), which is defined by (7ro(Xfc)0)(t) = tk4>{t) and {no(Yi)(f)){t) = 
^{t) for k,l = 1, . . . ,d. We will write W((i)+ for the set of all elements 
u G W{d)h such that {'Ko{u)(f),(j)) > for all G S{R'^). Clearly, 
^ W((i)^ C >V((i)+ while the converse is false by [301 Section 6]. 

For a given element c G yV{d)h write 



(5) 



inf c := sup{A G M | c - A ■ 1 G W{d)+} 
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if vro(c) is bounded from below. Otherwise write inf c = — oo. Clearly, 
inf c is equal to the infimum of the numerical range of 7ro(c), i.e. 

inf c = M{{no{c)v,v) \ v E S{R'^), \\v\\ = 1}. 

When 7Tq{c) is bounded from below, one can define the Friedrichs ex- 
tension 7ro(c)i7 of vro(c) and show that 

inf c = min a (tto (c) i?) , 

see e.g. [12] and the references therein. If cr(7ro(c)i?) is also discrete 
then inf c is equal to the lowest eigenvalue of 71q{c)f. 

Remark. There are well-known sufficient conditions on V{X) implying 
that the Schrodinger operator L = 7ro(— J2i=i + ^(^)) is bounded 
from below, is essentially selfadjoint (i.e. hp is the only selfadjoint 
extension of V) and has discrete spectrum, see [2DI Sections 8.1, 8.2]. 

Let c be an element of W((i)/i of even total degree degc. We propose 
the following method for computing inf c. Firstly, pick a sequence 
s = {sk)km of elements of W{d) with sq = 1. Secondly, solve the 
following sequence of "semidefinite programs" : 

/„x Find /ifc(c, s) = maxjr /i where J-" = {{fi, Z) | yU G M, 
Z e M„,(C), Z^ = Z, Z>0, bl{c-fi)bk = v^Zvk] 

where for every A; G N, 6^ = 11^=0 Vk is a vector of all monomials m in 
the generators of W{d) such that deg rri < ^ deg blcbk = deg &fc + 1 deg c 
and Hk is the size of Vk- 

Clearly, /ifc+i(c, s) > fik{c,s) for every k since bk divides bk+i in IQ. 
If 7ro(sfc) is invertible for every k, then, by (j5]) and ([6]), the sequence 
/ifc(c, s) is bounded above by inf c. The main question is what additional 
assumptions on c and s are needed to ensure that lim /ifc(c, s) = inf c. 
Our numerical experiments suggest that the only additional assumption 
required is that Sk are nonconstant but we are unable to prove that. 
What we can prove about convergence is summarized in Propositions 
2 and 3 below; see also Conjecture [H 

Recall that the leading symbol of an element 

deg c 

(in multiindex notation) is the element 

|a| + |;8|=degc 
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If c = c*, then c{X,i^) G R[X, For example, the leading symbol of 
N ■■= Eti -2 (X! - -l)e W{d)u isN = Eti -2 (Xf - ^f) . Note 
that N{X, tO = Eti 1 (Xf + ^f) is > if (X, ^ (0, 0). 

Proposition 2. Suppose that c G yV{d)h, 4| degc and c{X,i^) > for 
every X,^ eW^ with {X, ^) ^ (0, 0). Pick a G M+ \ N and a sequence 
rrik E'L in which every integer appears infinitely many times and write 
Sq = 1, Sk = N + {iTik + a) ■ 1 for k > 1. Then lim/ifc(c, s) = inf c. 

This is an immediate consequence of the fact that (t{N) = N and the 
following result of Schmiidgen which is a noncommutative analogue of 
Theorem [H (Assumption (1) can be relaxed slightly, see [HI Th. 1.2], 
but this does not help us here.) 

Theorem 2. [201 Th. 1.1] Suppose c G W{d)h, degc = 2m, satisfies 

(1) There exists e > such that c — £ ■ 1 G W((i)+. 

(2) c{X,i^) > for every X,^ G M*^ with (X,0 7^ (0,0). 
Finally, fix a E \ N and write M for the set of all finite products of 
elements N + {a + n) ■ 1, where n G Z. 

// m is even, then there exists b E J\f such that bob G J2 W('^)^- If^ 
is odd, then there exists b E JV such that Y^'^=i ^(^j + ~ ^ 

Proposition [3] is a variant of Proposition [2l We need some notation. 
We assume that d = 1 and write q = X and p = —lY. Every nonzero 
c G W(l) can be uniquely expressed as 

di d2 d2 di 

where 7^ and 7^ 0. In this case we say that c has multidegree 
{di,d2)- We fix two nonzero reals a and /3. Let iS be the monoid 
generated hy Si = p — ai, S2 = q — f3i and si, s^. It is an Ore set. 

Proposition 3. Suppose that c G W(l)h has multidegree (2mi,2m2), 
where mi,m2 G N, 72mi,2m2 7^ and f2m2 O'^'^d g2mi are positive on the 
real line. Since S is countable we can number its elements, say {ui)i^fi, 
assuming uq = 1. Write sq = 1 and let, for every k > 0, Sk+i G S 
be the common right multiple of Sk and Uk+i that exists by the Ore 
property. Then lim /ifc(c, s) = inf c. 

This is an immediate consequence of the following result of Schmiidgen: 

Theorem 3. [331 Theorem 5]. Let c be a nonzero hermitian element 
o/W(l) of multidegree (2mi,2m2), where mi,m2 G N. Suppose that: 
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(1) There exists a hounded self-adjoint operator T > on L^(M) 
such that 7ro(c) > T. 

(2) 72ml, 2m2 7^ and /2m2 '^'^'^ 5'2?tm orc positive on the real line. 
Then there exists an element s & S such that s*cs G ^ VV(1)^. 

Note that neither Proposition |2] nor Proposition |3] cover Schrodinger 
operators — + V{X) with polynomial potential V{X) of degree > 2. 
This case however fits into the following conjecture: 

Conjecture 1. Suppose that c G W(l)ft is bounded from below and 
Sk = aX + il for every A; G N. Then 2| degc and lim /ifc(c, s) = inf c. 

This conjecture is true if the following claim from (TU] is true: 

Claim. Let A be the algebra obtained from VV(1) with the addition of 
the generator {aX + il)'^ (for a G and the commutation relation 
[p, {aX + il)~^] = ia{aX + il)^^. Then every positive element c E A 
has a quadratic sum factorization c = ^ d^dk for some finite set of 
elements dk G A. 

The proof in [10] seems to have a gap (where they use a result of 
Schmiidgen). 

3. Implementation, Numerical Examples for d = 1 

The computation of fik{c,s) was implemented as follows. Firstly, 
the input for semidefinite programs (i.e. the matrices C, Ai and the 
numbers hi) was computed by Mathematica® (Wolfram Research) in 
rational (i.e. exact) arithmetics. Linear relations among Zij had to 
be solved before they were converted into the form tT{AiZ) = hi to 
ensure that the matrices Ai were linearly independent. Secondly, the 
input data was exported to Matlab® (Mathworks) where it was solved 
by either SeDuMi [55] or SDPT3 [M] semidefinite programming solver 
(through the Yalmip interface [U]) in floating point arithmetics. 

The problems of the basic method are illustrated by the following 
toy example in d = 1: 

Example 1. Write c = (2A^ + 1)^ = {X^-Y^f. Clearly, inf c = 1 but 
pretend we don't know that. The element c satisfies the assumptions 
of Proposition [2l We fix a = | and Sk = + a ■ 1. We will compute 
approximations /io(c, s), . . . , /i3(c, s) of inf c using sedumi and sdptS 
respectively. The results are in Table 1. The first approximation is 
very good but higher approximations are getting worse while the theory 
says they should be getting better. 
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k 


Hk{c,s) (sedumi) 


fik{c,s) (sdpt3) 


rik 







0.999999999760 


0.999999993360 


6 


15 


1 


1.000044013288 


0.999999989593 


15 


45 


2 


1.112894824977 


0.999999941936 


28 


91 


3 


73.69340728792 


24.89729311234 


45 


153 



Table 1. The table of Example 1. The numbers rik 
and ruk refer to the size and number of matrices Ai in 
semi definite programs used for computing //^(c, s). 

The lesson that we learn is that semidefinite programs should be 
kept as small as possible. The most natural way to do this is to exploit 
sparsity. (Later we will also discuss how to exploit symmetry.) For 
every 

c = 5^c„,;3X"F'^G W(l) 
write N'{c) for the convex hull of the set 

{{a — k, (3 — k) \ k = 0, . . . , min(a, /?)}. 

q,/3GN 

The point is that because of the relation YX — XY = 1 we must 
replace {(a, /?)} by {(a — fc, /3 — fc) | = 0, . . . , min(a, /?)}. As outlined 
in Section [H one can prove the property N'{fg) = N'{f) + N'{g) for 
every f,g & ^(1)- The property implies that the vectors Vk from 
the definition of fik{c,s) need not contain all monomials X"Y^ with 
a + (3 < degblfbk but only those from ^N'{bl{c — fi)bk)- This method 
works particularly well if bk depend only on X. 

Example 2. Let Ao(/3) be the lowest eigenvalue of vro(c/3) where 

= -Y^ + X^ + PX^. 

The values of Ao(/3) for various (3 were computed in [21 Table 1] to 15 
decimals. We refer to his values as "exact" . 

Let fik{(3) be the solution of the semidefinite program ([6]) for c = cp, 
bk = {i- X)'= and Wfc = (1, X, . . . , Y,XY,..., X^Yf. For each 

(3 e {0.0001, 1, 10000} and A; = 0, 1, ... , 14, we will compute by 
sdpt3. Finally, for each (3 we draw the semi-log plot of the sequence 
of relative errors of fik{/3) with respect to Ao(/3), i.e. the plot of the 
sequence k i-> log^g I '^"^^x^ip)^^^ I- The results are presented in Figure 
1. In theory these plots should decrease to — oo. By Figure 1, they 
decrease only during first 8-10 steps. Similar results are obtained for 
6fc = (1 + X2)'= and Vk = {l,X,..., X^''+^, Y, . . . , X^^Yf, see Figure 2. 
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Figure 1. The plots of k t-^ log^ 



I Ao(/3)-Mfc(/3)l 



for c 



'10 I Ao(/3) 

-Y'^ + + Sk = i — X and sparse f ^ (Example 2) 



yS^O.OOOl 

yS=l 

/3= 10000 



-4 - 




Figure 2. The plots of t-> log^o l^^^^^pf^l for c = 
-y2+x2+/3X^, Sfe = l + X^ and sparse Vk (Example 2) 



Whatever improvement we use, things eventually start to go wrong 
(because of rounding errors) and the question is when to stop. If we 
use sdptS, wrong values tend to undershoot, while if we use sedumi, 
they tend to overshoot. We can use this observation to formulate an 
empirical stopping criterium for sdpt3: If /io < /Ui < . . . < /x; > fii+i, 
then return as the best approximation for the lowest eigenvalue. We 
can also use min(/i; — — fii+i) as an estimate of its precision. 

There is no such stopping criterium for sedumi. We will use only sdptS 
in the sequel. 
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Remark. Another trick that sometimes improves stabihty in the com- 
mutative case is scahng X — AX, Y — > XY. In our case this does not 
work, because it violates the relation YX — XY = 1. On the other 
hand, the transformation X — )■ AX, Y X^^Y respects the relation 
but it does not improve stability. 

4. Exploiting symmetry - finite groups 

Suppose that G is a finite group acting on M.'^ by orthogonal trans- 
formations. This action induces in a natural way an action p on the 
polynomial ring M[X] and an action a on symmetric matrices that 
appear in our semidefinite programs, see [5]. The action a helps us 
put the matrices in our semidefinite programs into block diagonal form 
and thus reduce the amount of computation. It follows that every 
G-invariant sum of squares is a sum of squares of invariant and semi- 
invariant polynomials, see Theorem 5.3]. 

The same theory also works for Weyl algebras and even some more 
general algebras, such us enveloping algebras of Lie algebra. Finite 
generation of the ring of invariants for this situation was established 
in [1] by passing to the corresponding graded ring. An extension of O 
Theorem 5.3] to locally finite-dimensions actions of compact groups by 
*-automorphisms is provided by the last equation in the proof of [281 
Proposition 4]. 

Example 3. We would like to approximate the lowest eigenvalue of 

L = no{Y' + X') 

by exploiting symmetry. Let G = —1, —i} act on VV(1) by 

p(i)(X) = iY and p(i)(r) = iX. 

Clearly, p{L) = L. Since p{i){X + Y) =t{X + Y) and p(i)(X -Y) = 
—i{X — Y), it is more convenient to work with the generators 

X+Y , , X-Y 
a = ;= — and a = r= — . 

We start with the zero-th approximation, i. e. we would like to find 
the largest po such that Y'^ + X'^ — /iq ■ 1 is a sum of hermitian squares. 
We have to consider the 6-dimensional space M2 of all monomials of 
degree less or equal to 2. Eigenvectors of the restriction of p to M2 are: 

Ai = 1 : 1, a*a 
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The 6x6 matrices that appear in our semidefinite program can there- 
fore be assumed to be block diagonal with two 2x2 blocks and two 
1x1 blocks. We have therefore reduced the number of variables from 
15 to 8. Using sdptS, we get 

/io = 1.328427121 

To compute higher approximations, we need apropriate denominators 
bk such that h].Lhk is still G-invariant. Clearly, we can take for every 
G-semiinvariant polynomial, e.g. any element from A^. For 

hk = {2N + 1)'= = (X^ - = (2a*a + 1)^ 

we get 

/il = 1.396727721 

using the eigenvectors 

Ai = 1 : l,a*a,a^ (a*)^ (a*)V 

A2 = -1 : a^,{a*Y,a*a^,{a*)^a 

A3 = i : a, (a*)^, a*a^ 

A4 = — z : a*, a^, (a*)^a 

and 

/i2 = 1.396726593 

using the eigenvectors 

Ai = 1 : 1, a*a, a^, (a*)^, (a*)^a^, (a*)^a, (a*)^a^, a*a^ 

A2 = -1 : a^, (a*)2, a*a^, {a*ya, a^ {a*Ya\ (a*) V, (a*)^ 

X3 = i : a, (a*)^, a*a^, (a*)'*a, (a*)^a^, 

A4 = -2 : a*, a^, (a*)^a, 0*0"^, (a*)^a^, (a*)^ 

Our stopping criterium tells us that /ii is likely the best approximation 
we can get by this choice of bk- The method based on Conjecture [T] 
gives fi'^ = 1.396718666, /i'5 = 1.396726819, fi'^ = 1.396718409, i.e. a 
similar approximation and a similar estimate for precision. 

5. A CONJECTURE ABOUT RADIAL DIFFERENTIAL OPERATORS 

For every integer d > 1, we can identify the Hilbert space L^(R"''; r°'~^) 
with the subspace of L^(]R'^) consisting of radially invariant functions. 
Let iSrf(M"'") be the subspace of L^(R+; r'^~^) which corresponds to the 
space of radially invariant Schwartz functions on M"'. Let tt^ be the 
representation of the first Weyl algebra VV(1) which acts on iSrf(]R"*') by 

(7r,(X)/)(r) = r/(r) and (7r,(F)/)(r) = f{r). 

This representation is not a ^-representation however if we consider 
the standard involution on W(l) and the adjoint operation L 1— ?■ on 
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differental operators. To make it a *-representation we must consider a 
new involution L ^ L* := r'^-^L^r^~'^ which is conjugate to L ^ L"*". 

The aim of this section is to provide numerical support for the fol- 
lowing conjecture: 

Conjecture 2. If c E W(l)/i is such that iTd{c) > for some d then 
degy c is even and there exist A; G N and finitely many gi,hj G W(l) 
such that 

(1 + X'fc (1 + X' f = J2 9:9^ + E ^*3^hr 

i j 

Suppose we want to compute the lowest eigenvalue Aq of a dimensional 
radially invariant Schrodinger operator 



r'^ ^ dr dr' 

where V{r) is a real polynomial in r and 1/r. Pick the smallest m G N 
such that r'^'^~^d.-^v{r) has no negative powers. Conjecture [2] implies 
that the following sequence converges to Aq: 

jj,^ ■= sup{/i G M I 3gi, G W(l) : (1 + X^f {- X'^Y X^-^Y X"^ + 
+X2™+^-i(\/(X) - ^))(1 + X^)'^ = Y..glg, + h]Xh,}. 

We can write g*gi + = ul^Uk + vl{XB)vk, where Uk,Vk 

are suitable vectors of monomials and A, B are positive semidefinite 
complex hermitian matrices. Therefore, we can rewrite the definition 
of semidefinite program. 

Example 4. For = 1, . . . , 10 and 

V{r) = r + + r'^, 
we will compute approximations /io, • • • , of Aq using 



and 





^^2k+i+^^Y,...,X^^^'^YY rfodd, 
X^'^+i+i)^ deven. 



^ - ..,X2fc+f,F,...,X2'=-i+frf rfeven. 

and compare them with [6l Table 1]. The results are in Table 2. 

There is a problem with d = 1. Namely, by [6l Table 1], Aq is 
approximately 1.8306, while we can show (in exact aritmetics) that 
/io is approximately 1.9051. A possible explanation is that the eigen- 
function corresponding to 1.8306 is not in L^(M"'"). This has nothing 
to do with Conjecture [2J 
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d 


exact 


Un 








Ua 




2 


3.5644 


3.4973 


3.5623 


3.5643 


3.5644 


3.5643 


3.5630 


3 


5.3066 


5.2277 


5.3046 


5.3065 


5.3064 


5.3058 


5.3034 


4 


7.0746 


7.0073 


7.0730 


7.0746 


7.0746 


7.0677 


7.0405 


5 


8.8720 


8.8187 


8.8709 


8.8720 


8.8720 


8.8717 


8.8713 


6 


10.6987 


10.6488 


10.6978 


10.6986 


10.6986 


10.6874 


10.6030 


7 


12.5534 


12.5337 


12.5523 


12.5534 


12.5534 


12.5533 


12.5530 


8 


14.4348 


14.4229 


14.4345 


14.4348 


14.4342 


14.3901 


14.0928 


9 


16.3415 


16.3338 


16.3413 


16.3414 


16.3414 


16.3414 


16.3411 


10 


18.2720 


18.2664 


18.2719 


18.2720 


18.2712 


18.2041 


16.5905 



Table 2. The table of Example H For each d, the best 
approximation is underlined. Exact values for Aq are 
from [SI Table 1]. 



A 


exact 




Hi 


/U2 


f^3 


/U4 


At5 


0.0 


2.3381 


1.8899 


2.3193 


2.3368 


2.3380 


2.3300 


2.2971 


0.2 


2.1673 


1.7277 


2.1490 


2.1661 


2.1672 


2.1582 


2.1198 


0.4 


1.9885 


1.5583 


1.9706 


1.9874 


1.9884 


1.9783 


1.9372 


0.6 


1.8011 


1.3810 


1.7838 


1.8001 


1.8010 


1.7898 


1.7349 


0.8 


1.6044 


1.1951 


1.5878 


1.6035 


1.6043 


1.5915 


1.5447 


1.0 


1.3979 


1.0000 


1.3820 


1.3971 


1.3978 


1.3832 


1.3280 


1.2 


1.1808 


0.7949 


1.1657 


1.1801 


1.1807 


1.1641 


1.0800 


1.4 


0.9526 


0.5790 


0.9383 


0.9520 


0.9525 


0.9360 


0.8712 


1.6 


0.7127 


0.3516 


0.6992 


0.7121 


0.7124 


0.6937 


0.6241 


1.8 


0.4603 


0.1119 


0.4476 


0.4597 


0.4599 


0.4374 


0.3691 



Table 3. The table of Example O For each A, the 
best approximation is underlined. Exact values for Aq 
(rounded from 15 to 4 decimals) are from [Til Table I]. 



Example 5. Suppose that d = 1 and 

V(r) = h r. 

r 

We will compute approximations /io,---,/^5 of the lowest eigenvalue 
using Uk = {I,..., X^'^+i, Y,..., X'"'+^Yf and Vk = {1,..., X^^+Y 
and compare them with [TTl Table I]. The results are in Table 3. 

Example 6. Suppose that d = 1 and 

2 b 
V{r) = ar H — - 
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where a > 0, b > We apply one step of our method and divide the 
result by r on both sides. We get 



+ V{r)~2y^il + Jb + -)=g*g 



dr^ V 4 



where 9 = + y/ar ~ (| + \/b + j)r ^, which implies the inequality 



Ao>2v^(l + -^6 + ^). 
Similarly, if = 1 and 

V{r) = - + - 
where a < 0, b > —j, then we get as above 
d^ 



where = f - ^^/^^ - + ^/b + l) r-\ hence 



4 

An > 



il + 2^b + \y 

It is shown in [9j that both inequalities for Aq are in fact equalities but 
this is not clear from our method. If a > in the second case then 
Ao > because 



1 , /z, , 1 . 1 



d^ , a ,d ^ + ^^+%.^,d t + V^+i. 

+ V[r) = - + 



dr"^ r dr r dr r 

6. Other *-algebras 

The aim of this short section is to outline a possible extension of our 
theory from and W((i) to other *-algebras. 

Let A be a finitely generated real or complex unital *-algebra and 
IZ a family of equivalence classes of irreducible (possibly unbounded) 
^-representations of A. We can consider the elements of A as "poly- 
nomials" and elements of TZ as (evaluations in) "points" , see [32j . For 
every element c G A such that c* = c we can define 

inf c := sup{/i G M | 7r(c — /i • 1) > for every vr G TZ}. 

Clearly, our method for computing inf c can be applied to A if: 

• the monomials in the generators are linearly independent and 

• it satisfies an analogue of Theorem [H 
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Examples of such algebras are: 

• algebras of matrix polynomials [1], 

• enveloping algebras of finite dimensional real Lie algebras [3T] . 

• algebras of trigonometric polynomials [T8] . 

• finitely generated free real algebras [7], [12] (also [5]). 

Acknowledgement. 1 would like to thank Igor Klep for doing a part 
of programming and Konrad Schmiidgen for his comments. 
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